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ABSTRACT. In this paper we investigate a fractional order logistic map 
and its discrete time dynamics. We show some basic properties of the 
fractional logistic map and numerically study its period-doubling route 
to chaos. 
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C3 The concept of fractional order differentiation and integration is nearly 

as old as calculus itself Q. Already in the 17:th century Leibnitz made 

1-1 remarks on the fractional derivative of order 1/2 fl4l| . Despite many fun- 

damental results and important definitions found more than 150 years ago 
by among others Euler and Riemann, the field of applications of fractional 
calculus has mainly recently started to draw some interest |2ll5ll6ll7ll8lll2|. 

^ Applications to physics involve viscoelastic systems and electro-magnetic 

£S) waves |6], but more recently the dynamics of fractional order dynamical 

systems, involving fractional mechanics and fractional oscillators has been 

t— i studied HJ |12| . Fractional differential equations have for example arisen 

in Quantum Mechanics as well as in fractional generalizations of Einstein 
field equations ||5j [TTJ. Also fractional generalizations of various famous 
dynamical systems such as the fractional Lorenz-system, the Rossler sys- 
tems |6l and Chen's system have been investigated. In this paper we 
shall investigate the dynamics of a fractional generalization of the classic 
logistic map IfTol . 

2. Dynamics 

The discrete dynamics of unimodal maps is well known and has been 
extensively studied Ifl6l. In particular population dynamics often consider 
maps like the Ricker-family |Tl6| : 

(1) R\A X ) = Axe _ ^ x , A > 1 p > 0, 
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or the Hassel-family (See GH): 

(2) H KP {x) = A > 1 > 0. 

However the perhaps most famous unimodal map is the logistic map: 

(3) Q x (x) = Xx(l-x), A > 0. 

The logistic map is defined on the half-line x 6 [0, oo[, but all of its interest- 
ing dynamics takes place on the bounded interval / = [0, 1] for < A < 4. 
If we consider discrete time evolution by the mapping Q\ : I — »■ I, letting 
the state be x n at time n then we have x n+ \ = Q\{x n ) at time n + 1. In the 
prescribed interval the parameter A is increased which results in a period- 
doubling route to chaos. For a thorough background on unimodal maps 
and chaos see ||T6| . 

3. Fractional Calculus 

There are many models of fractional diff-integration, however the more 
popular ones are the Riemann-Liouville-type operator and Caputo-type 
operator [10], whereas the latter one is favored in cases of lack of initial 
conditions. We shall use the Riemann-Liouville fractional operator in order 
to construct the fractional logistic map. The Riemann-Liouville fractional 
integral is defined as follows. 

Definition 1. If f(x) G C([a, b]) and a < x < b then 

!_ r m 

T(a)J a (x-ty 

where a e] — oo, oo[, is called the Riemann-Liouville fractional integral of order 
a. In the same fashion for a e]0, oo[ we let 

which satisfies 

(6) D° a+ f(x) = I° a+ f(x) = f(x), 

and is called the Riemann-Liouville fractional derivative of order a. 

The fractional derivative is a special-case of the fractional integral. A 
fractional derivative of order 1/2 is called a semi-derivative and a fractional 
integral of the same order is called semi-integral |10fl. Both the fractional 
derivative and fractional integral satisfies the important semi-group prop- 
erty: 

Theorem 2. For any f e C([a,b]) the Riemann-Liouville fractional integral 
satisfies 

(7) 4Vf+/ w = &m 

fora>0,l3> 0. 



(4) J aV(*) == I 7Z \L a dt 
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A proof may be found in [10J. For more information regarding fractional 
calculus and its applications see 



4. Fractional dynamics 

4.1. Basic properties of the fractional logistic map. Many investigations 
of fractional dynamical systems contain fractional derivatives in time which 
is due to their continuous time derivative (See |HJ|6HZj])- There exists, how- 
ever, fractional differential equations in spatial sense, take the fractional 
Schrodinger equation for example J5l - In our fractional generalization of 
the logistic map we shall use a spatial fractional derivative. If we let the 
Riemann-Liouville fractional integral (|5) operate on the logistic map (|3) we 
get: 



A / 2r 
v ' v ' r a + 2) V a + 2 



x 1+a 



which is valid for a £ ]0, oof (for more details on fractional differentiation 
see 0), thus 

Definition 3. For all x > 0, A > and a £ ]0, oo[: 

(9) W*)--=ru^(l--^)* 1+a , 

1 (a + I) \ a + 2/ 

will be called the fractional logistic map (FLM) of order a. 

It is left as an exercise to show that the ordinary logistic map (|3| follows 
as a special-case when a = 0. The Riemann-Liouville fractional operators 
have special properties for the parameter- value a = 1/2 and are in these 
cases called semi-integral and semi-derivative respectively |9Ml0j[l4|. We 
shall therefore give the definition: 

Definition 4. For x > and A > the a = 1/2 order fractional logistic map 
appears like: 

and will be called the semi-logistic map. 

Due to the special property of fractional diff-integrals we can obtain the 
following useful theorem: 

Theorem 5. For all a, n e]0, oo[: 

in 

(id ^w=Qrw. 



Proof. It follows directly from definition [3] and the semi-group property 
(Theorem [2}. □ 
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Despite the fact that the fractional logistic map is valid for a G [0, oo] it is 
worth mentioning that probably the most interesting dynamics takes place 
for a G [0,1] since theorem [5] allows for translation to any other domain by 
iV-times integration/ differentiation for any N G Z. We can also give the 
following theorem: 

Theorem 6. For a G [0, 1]: 

(12) Q^(0) = Q^(l + ^a) = 0, 

holds for all A > 0. 

Proof. (0) = follows directly from definition [3] Also we find that 



Q"(l + |a) = by inserting x = 1 + ^a: 

v ' ^ M 2 1 r(a + 2)V a + 2/V 2 / 

Note that A > is a part of the definition of the fractional logistic map 
(Definition |). □ 

In order to find fixed points of the fractional logistic map we need the 
following criterion to be fulfilled: 

(14) Qx(x)=x, 
which is an equation on the form: 

(15) - + = 0. 

2 2a 

This equation has non-fractional solutions for a G [1, 2, 3, ...] which corre- 
sponds to simple integer integrations of the logistic map. For the fractional 
values a general solution is left as an open problem. Below we numerically 
study the dynamics of the fractional logistic map. 

4.2. Period doubling and chaos. With the aid of Mathematica 5 and soft- 
ware developed in [4] we were able to numerically study the fractional lo- 
gistic map. As an example we plot the fractional logistic map in the domain 

x G [0, 1], a G [0, 1]: 
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Fig.l. a G [0,1], x e [0,1]. 
Note that x = 1, a = 1 is a local maximum in the x-direction since the 
fractional logistic map for x = 1 and a = is zero due to the property 
of fractional diff-integrals in theorem [5] As an example we perform three 
iterates and end up with the following plot: 




Fig.2. a G [0, 1], x G [0,1]. 

In figure 2 we see Q®(Q^(Q®(x))). It is also worth mentioning that the gen- 
eral appearance for any choice of a G [0, 1[ appears similar to the ordinary 
logistic map (which is the special case when a = 0). Indeed, if we fur- 
ther examine the special case a = 1/2, the semi-logistic map, then we see 
that numerically a period-doubling route to chaos appears as A is increased 
from 4.5 to 6.1: 




Fig.3. a = 1/2, x G [0,3]. 

The appearance of this bifurcation diagram is similar to that observed for 
the logistic map [16|. Thus we have hinted that the fractional logistic map 
exhibits a period-doubling route to chaos. 

5. Discussion and open problems 

In this paper we studied the discrete time evolution of a fractional gener- 
alization of the logistic map; the fractional logistic map. We showed some 
basic properties of the fractional logistic map and numerically hinted a 
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period-doubling route to chaos for the special case a = 1/2. A rigorous 
proof of this is left for further investigation. The fractional parameter a in 
the fractional logistic map brings on a whole new range of possibilities for 
traditional discrete dynamics. The possibility of showing that "fractional 
chaos" arises when letting the traditional bifurcation parameter be con- 
stant while varying the fractional parameter is left as an interesting open 
problem. Furthermore, the dynamics of the fractional logistic map in the 
complex plane is also left as an open issue for further investigation. The 
connection between fractional discrete time dynamics, like the fractional 
logistic map, and other fractional dynamical systems such as the fractional 
R ossler and fractional Lorenz systems are also open issues. 
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